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I.  INTRODUCTION 

In  the  past,  ballistic  trajectory  determination  for  manual  or  computer  predicted 
ordnance  delivery  from  an  aircraft  was  determined  through  measurement  of  separation 
forces  on  the  piece  of  ordnance  via  wind  tunnel  or  captive  carry  measurements.  The  advent 
of  panel  method  codes  using  linearized  potential  theory,  such  as  A502i,  or  its  full  potential 
version  TranAir,  allow  for  a  cheaper  and  safer  method  of  predicting  separation  forces. 
Furthermore,  A502i  allows  for  any  arbitrary  configuration  to  be  modelled  within  the 
limitations  of  the  number  of  panels  and  networks  allowed  and  excluding  transonic  flow. 

The  purpose  of  this  work  was  to  determine  the  separation  forces  on  a  GBU-24 
carried  by  an  F-14  on  stations  3  or  6  or  both.  It  was  also  the  purpose  of  this  work  to 
provide  an  analysis  of  the  code  itself  to  see  if  it  is  a  viable  tool  for  the  study  of  flow 
characteristics  over  arbitrary  wing  configurations  for  use  in  the  Naval  Postgraduate 
School's  (NPS)  Department  of  Aeronautics  and  Astronautics.  The  majority  of  the  work 
was  conducted  on  the  NPS  computer  systems.  The  Department  of  Aeronautics  and 
Astronautics  Silicon  Graphics  Incorporated  (SGI)  workstations  were  utilized  for  most  of 
the  input  files  as  well  as  the  execution  of  the  code.  Due  to  the  amount  of  disk  space 
required,  storage  of  the  output  files  took  place  on  the  NPS  Computer  Center's  Y-MP  EL98 
Cray  computer.  The  bulk  of  the  GBU-24  data  was  calculated  using  the  SGI  workstations 
at  the  Naval  Air  Warfare  Center  in  Warminster. 

The  scope  of  this  analysis  was  to  understand  the  capabilities  of  the  A502i  code. 
The  approach  was  to  validate  A502i  against  existing  data  and  linear  theory.  The  code  was 
run  for  three  different  geometries  under  assorted  Mach  and  AOA  conditions.  Comparisons 
were  made  for  each  of  the  geometries. 


II.  OVERVIEW  OF  THE  A502i  CODE 

The  A502i  code  is  used  to  computationally  analyze  inviscid  subsonic  or  supersonic 
flows  about  arbitrary  configurations.  The  code  differs  from  other  panel  methods  in  that  it 
is  a  higher  order  panel  method;  that  is,  the  singularity  strengths  are  not  constant  on  each 
panel.  A502i  solves  the  linearized  potential  flow  boundary-value  problem  at  subsonic  and 
supersonic  Mach  numbers. 

The  aerodynamic  solution  provides  surface  flow  properties  (flow  directions, 
pressures,  Mach  number),  configuration  forces  and  moments,  sectional  forces  and 
moments,  and  pressures.  Additionally,  A502i  calculates  flow  properties  in  the  flow-field 
points  and  flow-field  streamlines.  Results  are  limited  to  subsonic  and  supersonic  cases 
(transonic  cases  excluded)  with  attached  flow.  Results  are  not  usually  applicable  to  cases 
where  viscous  effects  and  separation  are  dominant. 

A.  THEORY 

The  basic  equations  describing  the  flow  of  a  viscous,  compressible,  heat- 
conducting  fluid  are  the  Navier-Stokes  equations.  These  are: 
(a)  The  continuity  equation, 


dp  -        dp      i,   a(pVj) 

at         K  '    at    ~    ax, 


d         d         d 
where  V  =  ( — — ,  — — ,  — — )  is  the  gradient  operator  with  respect  to  the  location  vector 

x  =  (x,,x,,x3),  and  where  the  conventional  index  notation  is  used  instead  of  x=(x,y,z).    In 
addition,  t  is  time,    p(x,t)  is  the  density  and  V(x,t)  is  the  velocity  vector,  with  components 

V  =  (V,,V2,V3). 


( b)  The  conservation  of  momentum  equation, 

^(pV.Hl^tpV.V.l^+I^  +  pf,  0-1,2,3)  (2.2) 

where  x^  is  the  deviatoric  portion  of  the  stress  tensor  which  vanishes  for  a  frictionless 

fluid,  F(x,t)  is  an  external  body  force  per  unit  mass  exerted  on  the  fluid,  and  p(x,t)  is  the 
pressure. 

(c)  The  conservation  of  energy  equation 

^(pe+jplVf  +  pHI^KiW+iplVp+BJV,]  (23) 

where  e(x,t)  is  the  internal  energy  of  the  fluid,  k  is  the  coefficient  of  heat  conductivity  for 
the  fluid,  and  T(x,t)  is  the  temperature. 

(d)  The  equation  of  state 

f(p,p,T)=0  (2.4) 

where  the  function  f  depends  on  the  type  of  fluid.  The  derivations  of  these  equations  can 
be  found  in  Refs.  1,2  and  3. 

The  Navier-Stokes  equations  can  be  simplified  by  the  neglect  of  viscosity,  which  is 
equivalent  to  setting  the  deviatoric  stress  tensor  to  zero.  Combining  the  momentum  and 
continuity  equations  yields 

dVj         dp 
where  the  convective  derivative  operator  is  defined  as 

dt    at   f  'ax, 

Equation  (2.5)  is  referred  to  as  Euler's  equation.  The  continuity  and  energy  equations  can 


Equation  (2.5)  is  referred  to  as  Euler's  equation.  The  continuity  and  energy  equations  can 
be  reduced  to 

P  dt~( yl^r^-V-Vp+pV*?  (2-6) 

and  the  rate  of  increase  of  heat  per  unit  mass  is  given  by 

1  dp         d      1 

q  =  -V.(kVT)  =  -f +  p  — (-)  (2.7) 

p  dt         dt     p 

Equations  2.5,  2.6  and  2.7  can  be  reduced  to  a  single  equation  when  four  further 

assumptions  are  made.  First,  assume  isentropic  flow,  thus 

q=0  (2.8) 

Second,  assume  irrotationality 

VxV  =  0  (2.9) 

which  allows  for  the  introduction  of  the  potential  function  [Refs  2,3] 

V<I>=V  (2.10) 

Third,  assume  the  existence  of  a  freestream  potential  $,,  whose  gradient  is  the  uniform 

velocity  V^  attained  at  points  sufficiently  distant  from  the  disturbance  being  analyzed,  and 

thus  write 

(J>  =  4>-<I>  (2.11) 

and 

V  =  (u,v,w)  =  V4>  =  V<I)^+V(t)  =  Vjo+V(()  (2.12) 

The  quantities  <j>  and  v  are  called  the  perturbation  potential  and  velocity  IRef  3].  Fourth, 

assume  that 

|tf|2  «  a.  (2.13) 

everywhere,  where  a    is  the  freestream  speed  of  sound. 


Based  on  these  four  assumptions,  the  unsteady  potential  equation  is  obtained 
IRefs  1.3]: 

( 1  -  M2J(j>xx  +  <|>yy  +  $a  -  2M^xt  -  M>tt 

=  M'[-(y-  l)(2u+2<j)t+|v|2)V2(j) 

+  (2u  -  u2  )(|>xx  +  v:(t>yy  +  2vw(j)yz  +  w  2tya 
+  2(1  +  u)(v<|>xy  + w<))xz)  +  2(uut+ vvt+ wwt)] 

(2.14) 
Assuming  the  flow  conditions  do  not  change  with  time  yields  the  steady  non-linear 
potential  equation. 

=  M:[y(Y-l)(2u+|v|2)V2(f) 

+  (2u  +  u2  )<j)„  +  v2  <|)yy  +  2vw<J)yz  +  w 2  $a 
+  2(l  +  u)(v<j)xy  +  w<j)xz)] 

(2.15) 
When  Mx=0,  equation  (2.15)  reduces  to  Laplace's  equation, 

V2(j)  =  0  (2.16) 

For  the  case  of  Mx  *  0,  the  following  is  supposed, 

M2|v|«l-M2  (2.17) 

M2|v|«l  (2.18) 

which  are  small  perturbation  assumptions  [Refs.  1,2].  With  these  assumptions  the  steady 
non-linear  potential  equation  reduces  to  the  Prandtl-Glauert  equation  [Ref  1]: 

(l-M>„  +  <fryy  +  4>H-0  (2.19) 


Through  a  coordinate  transformation  [Refs.  1,2,3  j,  the  Prandtl-Glauert  equation  can  be  re- 
written as: 

s4>„  +  4>w  +  4>a-0  (2.20) 

where  when  s=l,  it  is  the  subsonic  case  and  Laplace's  equation  applies  and  when  s=-l,  it 
is  the  supersonic  case  and  the  wave  equation  applies.  Applying  Green's  third  identity  [Ref. 
1  ]  yields  the  following  integral  equation, 

<KP)  =  -—  m-if-lifi-V-i-ldS  (2.21) 

where  a  represents  the  source  strength  and  [i  represents  the  doublet  strength.    When 

supplemented  with  boundary  conditions,  it  is  equation  (2.21)  that  A502i  solves. 

A502i  solves  equation  (2.21)  through  a  discretization  process.  The  general  idea  of 
the  process  falls  into  two  parts.  The  first  is  developing  finite  dimensional  approximate 
representation  formulas  for  the  singularity  functions,  which  creates  a  system  of  equations 

with  unknown  coefficients,  Xt.  The  second  part  involves  solving  the  set  of  equations  for 
all  X,.  This  allows  for  completely  determining  the  source  and  doublet  functions.  Then,  by 

virtue  of  equation  (2.21),  the  potential  function  <|)(P)  is  determined  for  all  points  P,  solving 

the  problem. 

The  features  of  A502i  which  distinguish  it  from  predecessors  are  three-fold.  The 
first  is  a  feature  known  as  "continuous  geometry",  the  second  is  linear  source  and  quadratic 
doublet  variation,  the  third  is  continuity  of  doublet  strength. 

Most  panel  methods  approximate  the  configuration  geometry  with  panels  whose 
planform  is  a  quadrilateral.  Thus,  if  the  panels  themselves  are  planar,  only  a  small  class  of 
configurations  (such  as  cylinders  and  flat  wings)  can  be  described  without  gaps  being  left 
between  panels.  These  gaps  are  generally  small,  except  for  highly  twisted  surfaces.  The 


gaps  cause  little  numerical  error  in  subsonic  flow,  but  in  supersonic  flow,  the  cumulative 
effect  of  the  gaps  is  serious  [Ref.  1 1.  The  problem  is  not  associated  with  leakage  of  flow 
through  the  gaps,  but  with  the  doublet  strength  jumping  abruptly  from  a  non-zero  value  to 
zero  at  a  panel  edge  which  does  not  exactly  meet  the  adjacent  edge.  In  A502i,  gaps  are 
closed  by  means  of  panels  which  are  comprised  of  several  planar  regions. 

The  feature  of  linear  source  and  quadratic  doublet  variation  is  what  makes  A502i  a 
higher  order  panel  method.  The  basis  function  corresponding  to  a  source  parameter  is 
locally  linear,  while  the  basis  function  corresponding  to  a  doublet  parameter  is  locally 
quadratic.  This  is  what  allows  for  A502i  to  find  supersonic  solutions.  Numerical  solution 
of  the  wave  equation  is  far  more  sensitive  to  the  numerical  idiosyncracies  of  a  panel  method 
than  is  the  solution  of  Laplace's  equation  for  subsonic  flow.  Experimental  evidence  [Ref. 
1]  indicates  that  exact  surface  analysis  is  not  feasible  in  supersonic  flow  without  doublet 
continuity,  thus  the  potential  for  numerical  error  is  greatly  reduced  by  requiring  the  doublet 
singularity  strength  to  be  continuous  across  panels. 

B.  GENERAL  A502i  USAGE 

The  use  of  the  A502i  code  consists  of  generating  an  input  file,  which  can  be 
arbitrarily  named,  and  which  contains  the  information  defining  the  geometry  of  the 
configuration,  flow  field  points  of  interest,  the  flow  conditions  and  wakes.  The  process  of 
building  a  geometry  is  difficult  in  that  A502i  is  particular  about  its  input  format.  Simple 
configurations,  such  as  a  rectangular,  planform  wing  can  be  modelled  manually,  but  more 
complex  structures  require  a  pre-processing  program,  such  as  MACGS,  where  a  geometry 
can  be  graphically  built.  MACGS  will  output  a  data  file  in  a  format  that,  with  minor 
modifications,  via  another  pre-processing  program  that  can  move  the  data  from  three 
columns  to  six  columns,  will  be  readily  usable  by  A502i.  Currently,  the  school  does  not 
have  a  copy  of  MACGS,  but  it  can  be  acquired  through  McDonnell-Douglas.  To  complete 
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this  thesis,  MACGS  was  used  on  the  SGI  workstations  at  NAWC  Warminster.  Wakes 
also  must  be  constructed  in  the  same  manner  as  the  structure  to  be  analyzed.  More  detailed 
instructions  on  the  specifics  of  wakes  and  surface  geometries  can  be  found  in  Ref.  [4]. 
Appendix  A  is  a  portion  of  an  output  file,  but  lines  1  thru  1 1 20  are  an  exact  duplicate  of  the 
input  file. 

1.  Running  A502i  with  an  Existing  Executable 

Assuming  an  A502i  executable  file  (e.g.,  A502)  has  already  been  placed  in  a  user 
executable  directory  (e.g.,  /usr/local/bin),  the  only  other  necessary  items  needed  to  produce 
a  set  of  A502i  output  files  is  the  input  file  and  a  large  amount  of  storage  space.  Anything 
modelled  with  more  than  one  thousand  panels  total  will  use  more  than  one  hundred  mb  of 
disk  space.  If  the  maximum  number  of  panels  (20,000)  is  used,  the  disk  space  required 
will  be  on  the  order  of  2  gb. 

To  run  A502i,  enter  after  the  UNIX  prompt: 

A502  <input  file  >  output  file 
Prior  to  running  the  code,  it  is  highly  recommended  that  a  Cray  account  be  opened  and 
linked  to  the  department's  SGI  workstations.  This  is  done  by  assigning  the  same  user  i.d. 
number  to  the  Cray  account  as  is  assigned  to  the  account  with  the  department.  User  i.d. 
numbers  can  be  changed  by  the  computer  center  at  the  user's  request.  This  is  required  due 
to  the  limited  disk  space  available  to  individual  accounts  in  the  department.  Once  an 
account  is  opened,  log  on  to  a  department  workstation,  change  directories  to  an  existing 
Cray  directory,  for  example  (after  the  UNIX  prompt): 

cd  /jedi/dl /male tour 
Transfer  the  input  file  to  the  Cray  directory  and  execute  the  code.  The  screen  will  display 
what  portion  of  the  code  it  is  performing  and  how  long  it  took  to  perform  each  portion  in 
CPU  time.  The  code  outputs  numerous  files  in  addition  to  the  arbitrarily  named  output  file. 


The  two  output  files  of  interest,  in  the  vast  majority  of  cases,  are  the  arbitrarily  named 
output  file  and  the  ft  13  file.  In  order  to  run  another  solution  all  output  files  must  be  deleted 
or  renamed  prior  to  re-executing  the  code.  Relevant  results  should  not  be  kept  on  the  Cray 
account  as  files  on  disks  dl,  d2  and  ul  are  considered  temporary  storage  and  subject  to 
erasure  after  a  period  of  time. 

2.  Creating  an  A502i  Input  File 

The  input  file,  which  can  be  arbitrarily  named,  consists  of  two  portions,  the  largest 
being  the  geometry  data.  Appendix  A  is  a  complete  recreation  of  the  input  file  for  the 
GBU-24  with  canards.  The  file  begins  with  line  1,  $TITLE,  and  ends  with  line  1120, 
SEND.  The  line  numbers  are  for  reference  only  and  are  not  part  of  the  actual  input  file. 
The  first  portion  consists  of  creating  the  initial  conditions,  i.e.,  the  free-stream  Mach 
number  and  angle  of  attack,  the  type  of  analysis  to  be  performed,  i.e.,  solution  or 
datacheck,  what  types  of  output  that  are  to  be  included  in  the  arbitrarily  named  output  file, 
and  reference  points  to  be  used  in  calculating  forces  and  moments.  The  geometry  data 
consists  of  the  points  that  bound  each  panel,  that  in  turn  belong  to  a  specific  group  of 
panels  that  make  up  a  network.  The  overall  structure  being  modelled  consists  of  a  series  of 
networks.  A502i  can  run  up  to  150  networks  and  or  20,000  panels  with  a  limit  of  8,000 
panels  per  network.  Referencing  Appendix  A,  line  28  represents  the  first  network  of  the 
model,  a  canard.  Line  29  represents  the  number  of  networks  that  will  be  classified  under 
this  $POINTS  statement.  Line  30  indicates  what  kind  of  surface  the  network  will  be,  a 
three-dimensional  surface  with  flow  properties  to  be  calculated,  a  wake  and  a  base  are 
several  examples.  Line  31  is  the  number  of  y  points  and  the  number  of  x  points 
respectively  that  make  up  the  grid  of  that  network.  Line  32  is  where  the  panel  points  start. 
Reference  4  contains  detailed  instructions  on  the  options  and  meaning  of  each  of  the  non- 
geometry  inputs,  including  some  capabilities  not  shown  in  Appendix  A. 
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Two  types  of  solutions  can  be  run,  a  datacheck  and  a  full  solution.  Reference  4 
explains  how  to  enter  either  one  into  the  input  file.  The  datacheck  only  analyzes  the 
geometry.  This  can  be  accomplished  in  a  matter  of  seconds  for  a  simple  geometry  as  it  is 
only  running  the  first  several  portions  of  the  code.  The  full  solution  can  take  a  couple  of 
hours  for  a  geometry  of  the  size  of  4,000  panels.  The  datacheck  should  be  run  once  the 
geometry  has  been  modelled.  It  will  check  for  any  panel  edges  that  do  not  abut  properly, 
and  when  column  4  of  line  20  in  Appendix  A  is  a  1 ,  the  datacheck  will  list  the  unit  normal 
vectors,  which  must  be  facing  outward.  The  datacheck  will  also  see  if  the  wakes  are 
attached  properly.  A502i  is  capable  of  giving  warnings  both  on-screen  and  in  the 
arbitrarily  named  output  file  when  an  edge  or  a  wake  is  not  modelled  properly,  but  it  only 
lists  the  unit  normal  vectors.  The  directions  of  the  vectors  must  be  manually  checked  by 
the  user.  The  full  solution  performs  the  datacheck  first,  so  the  data  is  repeated  in  the 
arbitrarily  named  output  file.  Appendix  B  is  a  portion  of  the  output  file  that  contains  the 
summary  of  facing  surfaces.  Each  panel  edge  is  looked  at  to  see  what  it  abuts  against. 
Sections  such  as  wingtips,  leading  edges  of  a  flat  plate  or  any  surface  that  does  not  need  a 
wake  attached,  but  is  unabutted  to  any  other  panel  on  that  edge  will  draw  probable  error 
messages  or  warnings  from  the  code.  The  user  must  ensure  that  the  edge  is  not  supposed 
to  abut  against  anything  or  need  a  wake  attached.  If  that  is  the  case,  the  warnings  may  be 
ignored.  Appendix  D  is  the  first  page  of  the  portion  of  the  output  file  that  lists  the  unit 
normal  vectors.  The  three  columns  under  zc  are  the  x-y-z  coordinates  of  the  given  panel's 
center.  The  three  columns  under  znc  are  the  x-y-z  coordinates  of  the  unit  normal  vector.  In 
most  cases,  when  the  y  coordinates  are  of  the  same  sign,  then  the  unit  normal  vector  is 
pointing  outward. 
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C.  GEOMETRY  MODELLING 

Five  geometries  needed  to  be  modelled,  each  of  increasing  complexity.  Modelling 
proved  to  be  the  most  difficult  task,  in  that  A502i  is  a  FORTRAN  code  and  is  very  format 
sensitive,  but  the  sheer  number  of  points  that  need  to  be  generated  can  take  a  lot  of  time  and 
the  order  those  points  are  listed  in  the  input  file  is  what  determines  whether  or  not  the  shape 
is  correctly  modelled.  Of  the  five  geometries  modelled,  none  were  done  completely 
manually.  A  spreadsheet  was  used  for  generating  the  parabolic  arc  airfoil  and  the  deltawing 
since  those  structures  can  be  constructed  out  of  one  network,  excluding  wingtips  and 
wakes,  and  the  surface  can  be  defined  by  a  mathematical  function.  The  bombs  and  the  F- 
14  required  the  use  of  MACGS  to  be  properly  modelled.  MACGS  is  indifferent  as  to  the 
order  that  geometries  are  built,  and  often  doesn't  require  many  coordinate  inputs  if  building 
a  model  on  top  of  an  existing  IGS  file.  The  output  file  from  MACGS  is  automatically 
formatted  and  the  points  placed  in  the  appropriate  order  for  A502i  to  understand. 
Although,  the  order  may  be  reversed  where  the  unit  outward  normal  vector  is  concerned. 
MACGS  has  the  ability  to  output  files  in  several  different  panel  method  code  input  formats, 
including  PMARC.  Reference  4  gives  detailed  instructions  on  how  to  properly  order 
points  to  build  a  group  of  networks  that  will  model  a  geometry.  A502i  uses  a  right-handed 
coordinate  system  that  is  similar  to  an  aircraft  body  axes.  When  put  in  terms  of  a  wing,  the 
x  axis  is  positive  from  leading  edge  to  trailing  edge.  The  z  axis  is  positive  up  and  the  y  axis 
is  positive  out  the  right  wing 

1.  The  Parabolic  Arc  Airfoil 

The  parabolic  arc  airfoil  is  the  simplest  of  all  the  geometries.  The  airfoil  has  a  chord 
of  five  and  a  span  of  ten.  The  maximum  thickness  is  .15.  The  model  consists  of 
approximately  600  panels,  including  the  wake  and  wingtips.  A  spreadsheet  was  used  to 
develop  the  geometry  portion  of  the  input  file.  Line  32  of  Appendix  A  demonstrates  the 


12 


format  that  the  spreadsheet  used.  Rows  consist  of  two  points,  with  coordinates  x  1 ,  y  1 ,  z  1 , 
x2,  y2,  z2  using  a  format  of  6F  10.0.  The  chord  was  divided  into  25  points  (x  coordinate) 
from  trailing  edge  to  leading  edge  and  then  another  25  points  from  leading  edge  to  trailing 
edge  (bottom  half).  The  span  was  divided  into  12  points  (y  coordinate)  from  left  to  right. 
Due  to  the  symmetry  of  a  rectangular  planform,  the  y  coordinate  was  constant  along  the  50 
x  coordinates  that  constituted  a  chordwise  cross-section.  To  attain  a  maximum  thickness  of 
.15  the  formula, 

Z  =  .3*(A  __L)  (2.22) 

c      c 

was  utilized  to  generate  the  values  of  the  z  coordinates.  The  wingtips  simply  connect  the  x 
coordinate  on  the  top  side  with  it's  symmetrical  counterpart  on  the  bottom  side.  Due  to  a 
trailing  edge  composed  of  a  straight  line,  the  wake  is  modelled  by  a  single  panel  that  spans 
the  trailing  edge  and  has  a  length  aft  of  100.  Figure  2. 1  shows  the  panel  distribution  across 
the  top  surface  of  the  parabolic  arc  airfoil,  where  the  thickness  is  represented  by  the  color 
scheme.  A  panel  and  a  point  are  numbered  to  show  how  they  were  entered  into  the  input 
file. 

2.  The  Deltawing 

The  deltawing  represented  a  step  up  in  complexity  over  the  parabolic  arc  airfoil. 
The  chordwise  cross-section  is  parabolic,  while  the  spanwise  cross-section  is  linear.  The 
procedure  for  building  the  geometry  on  a  spreadsheet  was  the  same  as  that  for  the  parabolic 
arc  airfoil,  only  the  chord  length  is  not  constant  along  the  span.  For  simplicity  in  design, 
the  number  of  panels  per  column  of  panels  is  constant  on  the  deltawing,  as  on  the  parabolic 
arc  airfoil.  This  means  an  increasing  panel  density  in  the  direction  of  the  wing  tip.  The 
wake  is  modelled  the  same  as  the  airfoil.  The  right  wingtip  ended  in  a  point,  so  no  extra 
panelling  was  needed  to  close  any  gaps.    The  symmetry  of  the  deltawing  allowed  for 
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Figure  2. 1  Parabolic  Arc  Panel  Distribution 
further  simplification  and  reduction  of  the  code's  run  time  by  only  modelling  from  the 
centerline  to  the  right  tip.  A502i  allows  the  user  to  stipulate  whether  there  is  symmetry  in 
the  x-z  plane  and  or  the  x-y  plane  (see  line  5  and  6  of  Appendix  A).  This  means  that  the 
gap  between  the  top  and  bottom  panels  at  the  center  line  does  not  need  to  be  bridged  as  in 
the  parabolic  arc  airfoil  (symmetry  could  have  also  been  used  in  the  airfoil's  case).  The 
chord  of  the  deltawing  has  a  length  of  90  and  the  semi-span  has  a  length  of  15.  The 
maximum  thickness  occurs  midway  along  the  centerline  and  is  .05.  The  model  consists  of 
880  panels.  Figure  2.2  shows  the  panel  distribution  along  the  top  surface  of  the  deltawing. 
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Thickness  is  represented  by  the  color  scheme. 


Chord 


Figure  2.2  Panel  Distribution  of  Deltawing 


3.   GBU-24 

Wind  tunnel  experiments  were  run  on  the  GBU-24  without  canards  attached,  so  it 
was  deemed  relevant  to  build  a  model  with  and  without  canards  as  a  comparison  of  the 
code's  performance.  The  model  of  the  GBU-24  was  too  complex  to  build  with  a 
spreadsheet,  so  MACGS  was  used.  The  bomb  was  modelled  at  NAWC  Warminster  by 
superimposing  a  group  of  networks  on  top  of  an  IGS  file  being  displayed  by  MACGS. 
The  complete  configuration  consists  of  approximately  1300  panels.  Figures  2.3  and  2.4 
are  displays  of  the  GBU-24  with  canards,  with  Figure  2.4  including  the  wakes.  Figure  2.5 
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is  included  to  show  how  the  GBU-24  model  was  assembled.    Each  different  color 
represents  a  network. 

Several  features  of  the  geometry  are  relevant  to  point  out.  Two  of  them  are 
modifications  made  to  the  geometry  that  differ  from  the  actual  dimensions  of  the  bomb. 
Dr.  Alex  Cenko  of  NAWC  Warminster  has  extensive  experience  with  modelling  stores  in 
A502i.  The  modifications  were  made  on  his  knowledge  of  how  to  get  the  most  accurate 
results  from  the  code  when  modelling  stores.  The  first  is  to  model  the  fins  and  canards  as 
flat  plates,  i.e..  no  thickness,  which  A502i  allows  you  to  do  through  a  single  numerical 
change  in  the  input  code  for  each  network  that  represents  a  flat  plate  (see  line  30  of 
Appendix  A).  The  fins  and  canards  are  extremely  thin  when  compared  to  the  rest  of  the 
bomb,  and  to  add  a  third  dimension  to  the  geometry  complicates  the  construction  of  the  fin 
or  canard  for  several  reasons.  The  leading  and  trailing  edges  must  be  sharp  and  the  surface 
the  fin  or  canard  attaches  to  would  have  to  be  modified  to  abut  properly  with  two  edges 
instead  of  one.  Experience  has  shown  that  the  simpler  version  yields  accurate  predictions. 
A502i  is  an  inviscid  code,  so  it  cannot  take  into  account  separation  effects  on  its  own.  The 
GBU-24  does  not  have  a  flat  base.  In  reality,  it  is  more  bullet  nosed  in  shape.  However, 
at  the  speeds  with  which  the  bomb  is  being  analyzed,  separation  does  occur  near  the  trailing 
edge  of  the  bomb.  Experience  has  shown  that  truncating  the  end  into  a  flat  base  and 
designating  it  a  separated  flow  region  through  an  appropriate  input  (see  line  702,  column  1, 
Appendix  A)  yields  better  results  than  attempting  to  model  the  bomb  to  exact  physical 
dimensions.  The  last  feature  to  point  out  are  the  wakes,  as  seen  in  Figure  2.4.  A502i  has  a 
limitation  in  that  the  wakes  must  be  modelled  by  the  user,  and  they  have  the  same  abutment 
requirements  as  physical  surfaces.  Regardless  of  angle  of  attack,  the  wakes  remain 
stationary  with  respect  to  the  body  to  which  they  are  attached.  At  higher  angles  of  attack, 
the  wakes  are  no  longer  close  to  paralleling  the  free-stream  velocity.  Remodelling  the 
wakes  is  nearly  an  impossible  task.   The  fin  and  base  wakes  would  not  be  too  difficult 
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Figure  2.3  GBU-24  Geometry  (Wakes  not  Shown) 
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Figure  2.4  GBU-24  Geometry  (Wakes  Shown) 
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Figure  2.5  GBU-24  Geometry  (Network  Modelling) 
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because  they  do  not  abut  against  any  physical  surfaces  except  for  the  surfaces  they  are 
trailing  from.  The  canard  wakes  must  abut  properly  against  the  bomb's  body  all  the  way  to 
the  base.  Modification  of  those  wakes  would  entail  modification  of  the  entire  body,  or  in  a 
possible  simplification,  letting  the  wakes  remain  attached  along  the  body  until  the  base  of 
the  bomb  and  then  shifting  them  relative  to  the  free-stream. 

4.  The  F-14 

The  F-14  geometry  was  modelled  in  the  same  fashion  as  the  GBU-24.  The 
geometry  consists  of  approximately  1500  panels.  While  that  may  seem  fairly  coarse  for 
such  a  complex  structure,  experience  shows  that  it  is  all  that  is  required  to  get  accurate 
predictions.  The  primary  area  of  interest  is  the  underside  of  the  fuselage  forward  and 
between  the  two  nacelles.  Higher  panel  density  on  the  top  half  is  not  required.  Figures  2.6 
and  2.7  display  the  F-14  geometry  without  and  with  wakes  shown.  Several  omissions  are 
made  to  the  model  as  having  a  trivial  effect  on  the  analysis  or  no  effect  at  all.  Phoenix  rails 
and  bomb  racks  are  not  modelled  along  with  the  chin  pod  because  they  are  deemed 
insignificant  to  achieve  reasonably  accurate  predictions  over  small  angles  of  attack. 
External  tanks  were  not  considered,  but  could  be  modelled  much  in  the  same  way  as  the 
bomb  and  inserted  into  the  input  file  to  see  what  effects  the  drop  tanks  have  on  separation 
forces.  The  vertical  tails  and  horizontal  stabilizers  were  deemed  irrelevant  to  the  prediction 
of  the  separation  forces  and  were  left  out.  This  reduces  the  number  of  panels  and 
networks,  which  also  reduces  the  amount  of  time  it  takes  to  run  a  solution. 

5.  Combination  Geometries 

The  F-14  and  GBU-24  were  modelled  separately,  but  were  combined  together  as 
shown  in  Figure  2.8.  The  first  step  to  accomplish  this  was  using  the  FORTRAN  code 
NAVSEP  which,  among  many  of  its  functions  will  translate  coordinates  to  relocate  items  in 
the  flow-field.  Once  accomplished,  the  GBU-24  file  was  pasted  into  the  F-14  input  file. 
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Figure  2.6  F-14  Geometry  (Wakes  not  Shown) 
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Figure  2.7  F-14  Geometry  (Wakes  Shown) 
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Figure  2.8  F- 14  and  GBU-24  on  Station  3  Combined  Geometries 
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D.  GRAPHICS  VISUALIZATION 

One  of  the  pre/post-processing  codes  that  came  with  A502i  is  called  RAID. 
Currently  an  executable  exists  in  the  department's  computer  system  that  can  be  accessed 
through  typing  after  the  UNIX  prompt, 

raid 
Raid  is  a  basic  graphics  program  that  can  read  A502i  and  TRANAIR  input  files  and  display 
geometries  and  flow  properties  from  solutions.  After  accessing  raid,  it  will  ask  what  type 
of  input  file  it  is  being  asked  to  visualize.  It  can  handle  five  other  modes  of  input  besides 
A502i/TRANAIR  [Ref.  10].  It  will  then  ask  for  the  name  of  the  input  file.  A  prompt  will 
follow  asking  about  object  definition  matrices  which  needs  to  be  answered  by, 

EACH 
The  next  prompt  will  ask  if  the  panels  are  going  to  be  shaded  by  a  Cp  value.  Cp  is  used 
generically  in  that  Cp  can  be  displayed  or  Mach  or  any  of  the  relevant  49  surface  flow 
properties  [Ref.  4].  This  is  only  used  after  a  solution  has  been  run,  data  has  been  extracted 
and  a  colorscale  for  contour  plots  has  been  determined.  If  only  the  geometry  is  to  be 
displayed,  then  hit  carriage  return  for  all  the  next  questions  until  a  pink  window  appears 
with  a  menu  in  the  lower  left  corner  of  the  window.  An  anomaly  of  the  program  is  that  if 
you  want  to  display  the  wakes,  then  wake  display  must  be  deselected.  To  view  the 
geometry  select  view  on  the  menu.  All  selections  with  the  mouse  in  RAID  are  made  with 
the  center  mouse  button.  A  new  window  appears  with  a  menu  bar  at  the  bottom  and  left 
and  the  geometry  in  the  center.  From  there,  rotation,  translation,  scaling,  axes,  reflections 
and  other  manipulation  of  the  geometry  is  possible.  Figures  2.3  through  2.8  are  examples 
of  geometries  displayed  on  RAID. 

When  presentation  of  solutions  (i.e.,  Cp  or  Mach  contours)  are  desired,  the  use  of 
another  post-processing  program  is  required  to  generate  a  colormap  file.  The  program  is 
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called  crebar.  An  executable  currently  exists  on  the  department's  system.  Type  in, 

crebar 
after  the  UNIX  prompt.  The  program  will  ask  straightforward  questions.  Number  values 
associated  with  colors  available  can  be  found  in  Ref.  10.  The  color  file  can  be  saved  under 
any  name,  but  must  lie  in  the  same  directory  as  the  input  and  solution  files.  The  first  line  of 
the  color  file  will  list  four  numbers.  The  first  number  is  the  number  of  colors  assigned  to 
the  colormap  (248  maximum).  Occasionally,  the  color  bar  displayed  in  RAID  when  using 
a  colormap  will  disappear  when  certain  menu  items  are  selected.  To  prevent  this,  change 
the  last  three  numbers  to  read  6,  1,-1.  Plotting  outputs  from  RAID  requires  saving  the  file 
in  a  format,  such  as  RGB,  that  a  printer  will  recognize.  It  is  possible  to  change  the  text 
color  and  background  color,  the  default  is  black,  to  avoid  excessive  use  of  black  ink  in  hard 
copies. 
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III.  DATA  EXTRACTION 

When  a  solution  has  run  to  completion,  there  are  two  files  of  interest,  the  arbitrarily 
named  output  file  and  the  ft  13  file.  The  arbitrarily  named  output  file  contains  results  for 
everything  that  A502i  solves  for.  The  ft  13  file  contains  only  the  49  surface  flow  properties 
on  each  panel.  Appendix  D  is  the  solution  portion  of  an  arbitrarily  named  output  file  for  the 
first  network. 

For  purposes  of  displaying  flow  properties  on  RAID,  it  is  necessary  to  utilize  the 
ft  13  file.  A  post-processing  code  called  RAIDCONV  is  used  to  extract  the  specific 
information.  To  access  RAIDCONV,  type 

raidconv 
after  the  UNIX  prompt.  An  executable  currently  exists  in  the  department's  system.  The 
ft  13  file  must  be  in  the  same  directory  as  RAIDCONV  is  accessed  in.  RAIDCONV  will 
prompt  the  user  for  which  kind  of  panel  method  is  being  used  (A502i  is  one  of  three 
choices).  The  next  prompt  will  ask  for  the  name  of  the  ft  13  file.  The  last  prompt  will  ask 
for  the  flow  property  that  is  to  be  extracted.  A  file  called  ftl3.cp  will  be  created.  It  can  be 
renamed  for  purposes  of  multiple  flow  properties  extraction.  Abbreviations  for  the  49  flow 
properties  can  be  found  at  the  bottom  of  page  1  of  Appendix  D.  The  two  primary  flow 
properties  are 

LMACHU  for  local  Mach  number 
CP2ND  for  second  order  pressure  coefficient 
CP2ND  is  the  default  setting  for  RAIDCONV.  Once  the  ftB.cp  file  is  created,  RAID  can 
be  used  as  previously  discussed  to  display  the  flow  properties.     An  anomaly  of 
RAIDCONV  is  that  it  does  not  recognize  kt=20  type  wakes,  used  where  wakes  from  a 
wing  abut  against  a  body  [Ref.  4].   To  assist  in  extracting  all  the  data,  the  kt=20  wakes 
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should  be  placed  at  the  end  of  the  input  file.    In  general,  it  is  good  practice  to  place  all 
wakes  at  the  end  of  the  input  file  when  using  A502L 

The  arbitrarily  named  output  file  duplicates  the  data  found  in  the  ft  13  file  and 
includes  moments  and  forces.  A502i  will  sum  up  the  moments  and  forces  on  each  network 
and  for  all  networks  so  far  [App.  D].  The  moments  are  computed  based  on  the  coordinates 
entered  into  the  input  file  [Ref.  4  and  App.  A]. 
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IV.  RESULTS  OF  A502i  COMPUTATIONS 
A.  PARABOLIC  ARC  AIRFOIL  DISCUSSION 

This  simple  geometry  was  analyzed  primarily  to  evaluate  A502i's  capabilities  by  a 
comparison  to  known  linear  theory.  To  this  end,  the  geometry  discussed  in  section  II-C 
and  shown  in  Figure  2. 1  was  run  by  A502i  at  a  Mach  of  0.3  and  a  Mach  of  1 .5  at  an  angle 
of  attack  of  zero.  Two  of  the  49  flow  properties  that  A502i  computes  IRef.  4]  for  each 
panel  are  linear  Cp  and  second  order  Cp,  given  by 

CPLIN  =  -2uc  (4.1) 

CP2ND  =  -2iic-[(  l-M=c2)uc2+vc2+wc2]  (4.2) 

Where  uc,  vc  and  wc  are  the  compressible  components  of  the  perturbation  velocity. 
Figure  4.1  plots  the  linear  theory,  A502i  linear  and  second  order  results  for  the  subsonic 
case,  while  Figure  4.2  represents  the  supersonic  solution. 

Linear  theory  for  parabolic  arc  airfoils  is  outlined  in  Refs.  2  and  3.  The  equation 
representing  the  subsonic  case  is  given  by: 

1-x 


Cp(x) *Tr        ,  *(l-(.5-x)*ln 

x  *  chord  *  ^  1  -  NT 


)    where  0<x<  1  (4.3) 


Equation  4.3  includes  a  Prandtl-Glauert  compressibility  correction.    The  equation 
representing  the  supersonic  case  is  given  by: 

26 
Cp(x)=    ,  where  0  <  x  <  chord 

VM!-!  (4.4) 

1  2x     , 

and  0  =  xmai  *  (  — : — t t  ) 

™*       chord      chord - 
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Figure  4.1  Cp  Comparison  of  a  Parabolic  Arc  Airfoil  at  Mach  =  0.3. 
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Figure  4.2  Cp  Comparison  of  a  Parabolic  Arc  Airfoil  at  Mach  =  1.5. 


30 


Figure  4.1  shows  very  good  agreement  with  the  linear  theory  curve.  There  is  a 
small  but  noticeable  difference  between  A502i's  linear  results  and  second  order  results, 
with  the  second  order  results  being  more  accurate,  as  expected.  The  maximum  difference 
between  the  linear  theory  curve  and  the  second  order  A502i  curve  amounts  to  a  value  of 
2.5%  right  at  mid-chord.  The  gap  between  the  two  curves  from  .25  chord  to  .75  chord 
were  the  result  of  thin  panel  density  in  that  region. 

Figure  4.2  shows  excellent  agreement  with  the  linear  theory  curve.  The  A502i 
values  of  Cp  for  the  linear  and  second  order  analysis  are  virtually  identical.  There  are  small 
deviations  from  linear  theory  near  the  leading  and  trailing  edges,  but  this  is  expected  due  to 
numerical  error  associated  with  the  discontinuity  A502i  would  encounter  right  on  the 
leading  or  trailing  edges. 

B.  DELTAWING  DISCUSSION 

Reference  9  provides  data  on  Mach  distribution,  using  approximated  linear  theory, 
over  a  deltawing  of  the  configuration  discussed  in  section  II-D.  This  simple  geometry 
provided  another  test  of  A502i's  capabilities.  Figures  4.3  and  4.5  show  the  A502i  results 
for  the  subsonic  and  supersonic  case,  while  Figures  4.4  and  4.6  reflect  the  results  from 
Ref.  9.  For  both  cases,  good  agreement  is  found  with  the  linear  theory,  with  A502i's 
subsonic  analysis  being  physically  more  accurate  than  the  approximate  linear  theory,  while 
A502i's  supersonic  analysis  is  not  as  physically  accurate. 

A  comparison  of  Figures  4.3  and  4.4  reveals  several  points  of  interest.  The  Mach 
contour  representing  the  free-stream  value  is  given  by  the  dashed  line.  All  lines  outside  the 
dashed  line  represent  areas  where  the  Mach  value  is  less  than  free-stream,  and  inside  the 
dashed  line  is  where  the  Mach  value  is  more  than  free-stream.  The  location  of  where  the 
free-stream  Mach  contour,  in  Figure  4.3,  intersects  the  centerline  agrees  very  well  with 
Figure  4.4.   However,  Figure  4.4  does  not  have  the  contour  extending  all  the  way  to  the 
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tip.  This  is  a  physical  limitation  of  the  approximate  theory  used  in  Figure  4.4  and  A502i  is 
giving  a  more  realistic  solution.  Figure  4.4  suggests  that  the  peak  Mach  value  occurs  at 
approximately  two-thirds  chord  along  the  centerline.  The  A502i  results  show  the  peak 
Mach  contour  occurring  out  midway  along  the  semi-span.  Those  Mach  values  are  less  than 
1%  larger  than  the  yellow  Mach  contour  surrounding  it,  and  can  be  attributed  to  how  the 
panel  density  increases  with  movement  towards  the  wingtip.  A502i  performed  very  well 
for  this  subsonic  case. 

A  comparison  of  Figures  4.5  and  4.6  shows  that  A502i  did  not  perform  as  well  as 
in  the  subsonic  case.  Again,  the  contour  representing  the  free-stream  value  of 
Mach  =1.414  is  given  by  the  dashed  line.  All  lines  forward  are  below  free-stream  and  all 
lines  aft  are  above  free-stream.  Figure  4.6  shows  the  intersection  of  the  free-stream  Mach 
contour  on  the  centerline  occurring  at  approximately  39%  chord,  which  is  in  excellent 
agreement  with  A502i's  result.  Figure  4.6  shows  the  peak  Mach  value  occurring  at  the 
trailing  edge  on  the  centerline.  This  makes  more  physical  sense  than  the  results 
that  A502i  yielded.  The  maximum  thickness  of  the  deltawing  occurs  along  the  centerline, 
allowing  for  greater  expansion.  The  discrepancy  may  be  attributable  to  panel  density  and 
accumulation  of  numerical  error.  A  close  study  of  the  A502i  results  reveals  some 
discontinuities  along  the  column  of  panels  out  at  the  wing-tip  which  would  have  adversely 
affected  the  solution  and  caused  errors  to  propagate  along  the  semi-span. 


32 


Lines  of  Constant  Mach  (Minf=01) 


Semi-Span 
Figure  4.3  AS02i  Mach  Contour  Plot  (M«,=.01) 
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Figure  4.4  Approximate  Linear  Theory  [Ref.  9]  Moo=0 
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Lines  of  Constant  Mach  (Minf=1 .414) 


Semi-Span 


Figure  4.5  A502i  Mach  Contour  Plot  (Moo=1.414) 


Figure  4.6  Approximate  Linear  Theory  [Ref.  9]  Moo=1.414 
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C.  GBU-24  FREE-STREAM  (NO  CANARDS)  DISCUSSION 

A  free-stream  measurement  of  separation  forces  on  the  GBU-24,  without  canards, 
was  conducted  in  a  wind  tunnel  for  various  Mach  numbers  from  .8  to  1.2  JRef.  7J.  Since 
A502i  uses  linear  potential  theory,  the  model  of  the  GBU-24,  without  canards,  was 
evaluated  at  both  Mach  .8  and  1.2,  avoiding  the  transonic  regime,  to  ascertain  the  accuracy 
of  the  code  with  the  given  geometry.  Normal  forces  and  pitching  moments  for  both  cases 
are  plotted  and  compared  to  the  wind  tunnel  data. 

1.  Subsonic  Case  (Mo>=0.8) 

The  GBU-24  model,  without  the  canards,  was  run  for  angles  of  attack  varying 
from  -10  to  +10  degrees  in  two  degree  increments.  Values  much  higher  than  that  ran  into 
wake  modelling  problems  as  the  wake's  angle  relative  to  the  free-stream  was  getting  large 
enough  that  results  would  become  questionable,  and  remodelling  the  wake  was  too  difficult 
for  such  a  complex  geometry.  The  results  of  the  A502i  analysis  are  displayed  in  Figures 
4.7  and  4.8.  For  angles  of  attack  between  -4  and  +4  degrees,  A502i  does  a  good  job  of 
predicting  the  separation  forces.  The  pitching  moment,  which  happens  to  be  unstable 
without  the  canards,  is  approximately  linear  over  the  -4  to  +4  degree  range  and  is  the 
limiting  factor  to  the  models  accuracy.  The  normal  force  is  approximately  linear  over  a 
wider  range,  and  A502i  does  a  good  job  of  predicting  the  normal  forces  from  -6  to  +6 
degrees. 

2.  Supersonic  Case  (Moo  =1.2) 

The  results  for  the  subsonic  case  demonstrated  that  the  effective  range  of  angle  of 
attack  that  A502i  needed  to  explore  was  from  -6  to  +6  degrees.  Figures  4.9  and  4. 10  plot 
the  comparison  of  wind  tunnel  data  versus  A502i  results  for  pitching  moment  and  normal 
force.  The  results  for  the  supersonic  case  are  slightly  better  than  that  of  the  subsonic  case. 
The  actual  pitching  moment  of  the  GBU-24  is  approximately  linear  over  a  wider  angle  of 
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Figure  4.7  Comparison  of  Pitching  Moments 
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Figure  4.8  Comparison  of  Normal  Forces 
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attack  region,  but  fluctuations  in  the  data  at  -6  and  +6  degrees  means  that  the  model  is  still 
only  viable  from  -4  to  +4  degrees.  The  normal  force  line  is  nearly  linear  from  -10  to  +10 
degrees  and  extrapolating  the  A502i  results  out  to  10  degrees  would  still  yield  good 
predictions. 
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Figure  4.9  Comparison  of  Pitching  Moments 
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GBU-24  Free-Stream  (M=1.2) 
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Figure  4.10  Comparison  of  Normal  Forces 

D.  GBU-24  FREE-STREAM  (WITH  CANARDS)  DISCUSSION 

As  in  the  case  with  no  canards,  free-stream  measurements  of  the  separation  forces 
on  GBU-24  were  taken  from  Mach  .8  to  1.2  in  a  wind  tunnel  |Ref.  7|.  Again,  due  to  the 
limitations  of  linear  theory  inherent  in  the  code,  an  analysis  was  done  for  Mach  numbers  of 
.8  and  1.2  to  minimize  transonic  effects.  Even  with  the  more  complex  geometry,  A502i 
does  an  accurate  job  of  predicting  the  separation  forces  over  the  range  of  angles  of  attack 
that  are  approximately  linear. 

1.  Subsonic  Case  (Moo=0.8) 

The  GBU-24  model,  with  canards,  was  run  in  two  degree  increments  of  angle  of 
attack  from  -10  to  +10.  The  wake  modelling  limitation,  as  well  as  reviewing  the  data  from 
the  wind  tunnel  measurements  |Ref.  1\  showed  the  non-linearity  of  the  separation  forces  at 
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the  higher  values  of  angle  of  attack,  precluded  any  attempts  to  predict  forces  beyond  the 
aforementioned  angle  of  attack  interval.  The  results  of  the  A502i  analysis  are  displayed  in 
Figures  4. 1 1  and  4. 12.  The  addition  of  the  canards  makes  the  pitching  moment  stable,  but 
linear  over  a  smaller  region  than  without  the  canards.  A502i  gave  accurate  results  from  -2 
to  +2  degrees  angle  of  attack  when  predicting  pitching  moment.  The  prediction  of  normal 
forces  fared  better,  showing  accurate  results  from  -3  to  +3  degrees  angle  of  attack. 
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Figure  4.11  Comparison  of  Pitching  Moments 
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Figure  4.12  Comparison  of  Normal  Forces 

2.  Supersonic  Case  (Moo=1.2) 

As  in  the  case  with  no  canards,  the  region  of  accuracy,  with  the  model  used,  was 
assumed  to  be  less  than  +  or  -  10  degrees  angle  of  attack.  Cases  were  run  from  -6  to  +6 
degrees  angle  of  attack  in  two  degree  increments.  A  comparison  of  A502i  results  with 
wind  tunnel  data  is  shown  in  Figures  4. 13  and  4.14.  For  both  the  pitching  moment  and  the 
normal  force,  A502i  does  a  much  better  job  of  prediction  then  when  subsonic.  The  wind 
tunnel  data  is  nearly  linear  in  both  pitch  moment  and  normal  force  from  -8  to  +8  degrees 
angle  of  attack.  Extrapolating  the  A502i  data  out  to  +  or  -  8  degrees  angle  of  attack,  shows 
excellent  agreement  with  the  wind  tunnel  data. 
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Figure  4.13  Comparison  of  Pitching  Moments 
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Figure  4.14  Comparison  of  Normal  Forces 
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E.  F-14  DISCUSSION 

The  GBU-24,  without  canards  was  located  at  station  3  of  the  F-14  via  the  NAVSEP 
code.  In  the  process  of  running  a  solution,  the  combined  geometries  were  found  to  have  a 
total  of  150  networks  (A502i's  maximum).  The  combined  geometry  had  about  4,000 
panels,  far  short  of  the  20,000  panel  maximum,  so  there  was  room  for  more  detailed 
modelling,  but  there  was  not  a  chance  to  insert  the  GBU-24,  with  canards,  into  the  F-14's 
flow  field  and  analyze  the  forces  on  the  bomb.  The  geometry  with  the  canardless  bomb  ran 
to  a  solution  that  appeared  to  be  valid,  so  there  is  a  high  degree  of  confidence  that  if  the 
number  of  networks  could  be  reduced  to  allow  the  GBU-24,  with  canards,  to  be  inserted 
into  the  F- 14  flow  field,  the  code  would  yield  accurate  predictions  at  small  angles  of  attack 
on  the  forward  stations.  To  reduce  the  total  number  of  networks  by  combining  existing 
networks  would  have  required  a  large  time  investment  and  the  use  of  MACGS,  which  the 
department  currently  does  not  possess,  the  two  reasons  why  it  was  not  done.  Figure  4.15 
shows  a  Mach  distribution  of  the  solution  of  the  canardless  bomb  and  F-14  at  Mach  =  0.8 
and  0  degrees  angle  of  attack 

F.  POST-PROCESSING  DISCUSSION 

The  Mach  values  for  the  subsonic  case  of  the  GBU-24,  with  canards,  at  4  degrees 
angle  of  attack,  were  extracted  from  the  ft  13  file.  These  values,  used  in  conjunction  with 
RAID  are  shown  in  Figures  4.16  and  4.17.  The  color  distribution  over  the  nose  in  Figure 
4.16  indicates  that  the  bomb  is  at  an  angle  of  attack,  and  scanning  the  rest  of  the  model 
showed  no  discontinuous  solutions,  which  is  generally  represented  in  A502i  by  a  Mach 
value  of  0  or  1,000.  The  visual  representation  is  a  quick  way  of  telling  if  A502i  ran  an 
accurate  solution.  The  only  other  way  is  to  individually  check  the  Mach  or  Cp  values  of 
each  panel  in  the  ft  13  file  or  the  arbitrarily  named  output  file.  The  other  point  of  interest  in 
Figure  4. 16  is  the  lack  of  panel  density  along  the  mid-section  of  the  bomb.  The  goal,  in  the 
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case  of  stores  separation  prediction,  is  to  have  as  simple  a  model  as  possible  that  still  gives 
accurate  predictions.  The  fewer  the  number  of  panels,  the  shorter  the  run  time.  The  fact 
that  A502i  is  a  higher  order  panel  method  allows  the  luxury  of  using  fewer  panels.  Figure 
4.17  highlights  the  approach  used  to  take  into  account  separation  effects  as  discussed  on 
page  16. 
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Figure  4.15  Mach  Distribution  over  GBU-24  and  F-14     (Mco=0.8,    alpha=0°) 
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Figure  4.16  Mach  Distribution  over  GBU-24  (M«>=0.8,    alpha=4°) 
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Figure  4.17  Mach  Distribution  over  GBU-24  (Moo=0.8,    alpha=4°) 
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V.  SUMMARY  AND  CONCLUSIONS 

The  main  goal  of  this  analysis  is  to  determine  the  accuracy  of  A502i  on  both  simple 
geometries  and  complex  geometries.  To  accomplish  this,  A502i  is  compared  with  results 
available  from  linear  theory  and  wind  tunnel  experiments.  This  allows  conclusions  to  be 
drawn  on  the  capabilities  as  well  as  the  limitations  of  A502i. 

In  general,  A502i  can  accurately  predict  flow  properties,  forces  and  moments  on 
simple  and  complex  geometries  at  low  angles  of  attack.  The  predictions  are  valid  over  a 
wide  Mach  range,  from  0  up  to  and  including  0.8  and  from  1.2  and  above.  The  supersonic 
solutions  are  available  due  to  A502i's  higher  order  capabilities. 

The  limitations  of  A502i  are  consistent  with  most  panel  methods.  A502i  cannot 
predict  flow  dominated  by  viscous,  separated  or  transonic  effects.  It  cannot  predict  flow 
with  different  total  pressures,  such  as  flow  properties  inside  a  jet  plume  or  a  propeller 
slipstream  swirl.  The  biggest  shortcoming  of  A502i  is  its  inability  to  handle  unsteady  cases 
and  automatically  determine  wake  shapes. 

Experience  or  knowledge  of  the  flow  properties  around  the  geometry  being  tested  is 
important  in  building  an  accurate  model.  An  accurate  model  may  not  be  physically 
accurate.  Right  test  results  revealed  a  yawing  moment  on  the  GBU-24  that  was  not 
discovered  in  wind  tunnel  experiments  when  the  bomb  was  carried  on  an  aft  station.  The 
yawing  moment  may  be  caused  by  the  fact  that  the  canards  are  not  fixed,  but  spring- 
damped.  A502i  predictions  would  not  be  accurate  without  inputting  a  moment  to  simulate 
the  deflection  of  the  canards,  since  the  canards  are  fixed  by  the  geometry.  A502i,  while  a 
powerful  tool  in  terms  of  cost  savings  and  time,  cannot  completely  substitute  for  wind 
tunnel  experiments  and  flight  tests,  as  constructing  a  complex  geometry  to  exact  physical 
specifications  will  probably  not  yield  accurate  predictions. 
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APPENDIX  A.  GBU-24  OUPUT  FILE  (INPUT  FILE  PORTION) 
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APPENDIX  B.  GBU-24  OUPUT  FILE  (EDGE  ANALYSIS) 
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